Superconducting correlations in any dimensionality. 
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We consider the static anomalous Cooper loop for the electron gas of arbitrary spatial dimen- 
sionality, D. This object enters the mean-field equation for the superconducting temperature, T c . 
The closed expression in r— space is found at T = as an analytical function of D. Its counterpart 
in q— space has a logarithmic singularity at q = and no singularities at q = 2&f for any D. We 
extend our analysis for the case of finite T and note that in two dimensions the bound state of two 
electrons is possible even in the limit &f — > 0. We comment on the possible implications of this fact 
for the theory of high-T c cuprates. 
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The superconductivity is one of the collective phenom- 
ena in the itinerant electron systems which attracts much 
of the theoretical interest. Last decade this interest was 
strongly motivated by the puzzles of the high-T c cuprates. 
The basic ingredient of these materials is the copper oxide 
planes, which let one think that the interesting physics 
of cuprates is connected with the reduced dimensionality 
of the system, D = 2. 

It is trivial nowadays to say that the superconductiv- 
ity corresponds to the formation of the bound states of 
electrons with opposite momenta (Cooper pairs). EJ One 
may also say here, that given an appropriate sign of the 
quasiparticle interaction <?, the onset of the superconduc- 
tivity is revealed as the singularity in the RPA-collected 
series of diagrams (Cooper ladder) as T approaches T c 
from above. This divergency means the instability of the 
initially assumed ground state and makes possible the 
non-zero superconducting (SC) order parameter. The 
transition temperature T c could be found from the equa- 
tion <?II(q = 0,T C ) = — 1. The properties of the static 
Cooper loop II (q) thus define T c on the mean- field level 
of consideration. 

Generally, there are two known types of instabilities 
in the itinerant electron systems. Similarly to the super- 
conductivity and the role of the Cooper loop for it, one 
also discusses the magnetic instability (spin-density-wave 
states) and the corresponding magnetic RKKY loop. The 
analysis of both types of instabilities ia.,Darticularly im- 
portant for the one-dimensional metalsaO. 

In the recent paperu we presented the method which 
enabled us to find in a simple manner the magnetic 
RKKY loop for the spherical Fermi surface (FS) and any 
dimensionality D of the system (see alscS). The closed 
expressions were found both in r— and q— representa- 
tions. The suggested method also helped us to obtain 
the closed analytic expressions for the RKKY loop in a 
metal with strongly anisotropic FS.El 

In the present paper we extend our analysis for the 
case of superconducting correlations in a metal of arbi- 
trary spatial dimensionality. We report the closed ex- 
pression for the Cooper loop in the limit of zero temper- 
ature for any D in r— representation and its counterpart 



in q— space. In the two-dimensional (2D) situation the 
Fermi momentum disappears from the final expressions. 
It means that the tendency to superconducting pairing 
remains in the limit of vanishingly small 2D Fermi surface 
; it is fairly distinct from the case of three spatial dimen- 
sions. It is interesting to note that the Fermi surface in 
the form of small 2D pockets was recently discussed in 
connection with the high-T c cuprates.Q Thus the present 
study might be useful for the theory of the high-T c com- 
pounds, and we briefly comment on it. 

We begin by defining the basic object of our consider- 
ation. The static Cooper loop in r— representation could 
be written as 



n(r)=T^GKr)G(-M,r), 



(1) 



with Matsubara frequency u>i = irT(2l + 1) and the elec- 
tronic Green's function given by 



G(iw,r) = 



d D \t exp(ikr) 

(2tt) d iu> — £k 



(2) 



We focus our attention below first at the case of low 
temperatures when one can use the limiting relation 
T^2 t — > J dio / (2n) . We also explicitly assume the 
quadratic electron dispersion in D dimensions : 

e k = k 2 /2m- fi, (3) 

with the Fermi energy p — k 2 F /2m. Using the definitions 

v = D/2 — 1, z = fj, + iuj, p = 2mr 2 , 

the Green's, function could be conveniently represented 
in the forma : 



G(fc,,r) = -(f 



K^^Tp). (4) 



In this equation the branch of the root \J—zp should 
be chosen from the condition of its positive real part. In 
particular, this latter condition means that the argument 



1 



of Mcdonald function K v (y/—zp) has a discontinuity at 

LO = 0. 

The R— dependence of (jlLcould be evaluated with the 
use of Macdonald's formulaS, namely 
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(5) 



In the considered situation we have a, b = 
y/—p(lJ, ± iuj) and the product ab in (||) is continuous 
and real at u> = 0, although each of the arguments a, 6 is 
not. In particular, it guarantees the applicability of Eq. 
(|), defined by the condition! |Arg(a + b)\ < tt/4. Us- 
ing the integral representation of the appearing function 

K v P\/ H 2 + oj 2 /t and performing first the Gaussian in- 
tegration over uj, we obtain after simple calculation 

n(r) = [iN 2 F T 2 [v + 1] <f>(k F r), 

with the density of states at the Fermi level 



N F = 



27rl> + l] V 47r 
The range function entering (fl) is found in the form 
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ds cosh(2i/ + 1)5 e 



-2x sinh i 



2[2/x] 2lJ+l S , 2v+1 (2x), 



(0) 
(7) 

(8a) 
(8b) 



with the Lommel's function S a .b(x). It follows from (||) 
that at large distances II(r) oc r~ D and at small dis- 
tances II(r) oc r 2 ~ 2D . In the even dimensions $(x) is 
expressed through the Neumann polynomials Op>-i(2x); 
for the odd dimensions $(x) is reduced to Bessel func- 
tions. It is interesting also to note that (||) implies that 
IT(r) satisfies the following differential equation : 

d 2 U(r) 2D - 1 dU(r) ,, , . /fe\ D ,„ v 

— + -^+4fc|nr =4m ( (9) 

ar z r ar \ Zirr J 

In the particular cases D = 1,2,3 we have!: 



n(r) = -[H (2fc i rr)-y (2fc i rr)], D = 1, 



(27rr) 2 ' 



D = 2, 



(10) 
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H 2 (2fc F r)-r 2 (2fc i ,r) 
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Note that in the two-dimensional case the value of k F 
does not enter n(r) and Np. 

Knowing the expression (|^) for the correlations in 
R— space one finds its correspondence in q— space as fol- 
lows. 



n(<j) = J d D re lclr U{r) 



2ixr 



is+l 
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A straightforward application (cf. below) of the relevant 
formulasa shows that II (q) is expressed via the Gauss 
hypergeometric function 2-F1 [a, b, c; z]. Introducing the 
dimensionless momentum n — q/(2k F ) we have 



n(«) 



Np 



2vk 2 
Np 
"2V 



F[1,1;1-j/;1-k- 2 ] 



F 1, —v; 1 — v;l — k 



(12a) 
(12b) 



where ([l2"|b) analytically continuates (|l2|a) onto the re- 
gion hk > 1. It means that, in contrast to the RKKY 
loopp the Cooper loop H(q) reveals no singularity at 
q = 2k F , which is verified also by the non-oscillatory 
character of II (r) in (||). The logarithmic divergence at 
q = in (y_2|) corresponds to the aforementioned decrease 
n(r) ~ r at large distances. 

Actually, the situation with ( [l2] ) is more delicate. We 
managed to obtain the expression analytically depen- 
dent on D, despite the fact that at D > 3 one has 



n(r) ~ r 2 ~ 2D a t small r and (|il| ) diverges at the lower 
limit. While in even dimensions one has v integer and 
( |l2| ) divergent, in odd dimensions 11(g) has a finite value. 
To resolve this paradox, we let v be complex and use 
the intermediate representation of (jll|) in the form of the 
Mcllin-Barnes integral. □ It could be seen then that the 
result ( [l2] ) corresponds to the subtraction of the first few 
divergent terms proportional to r 2 ~ 2D ', r~ 2D , . . . , r ~ x ~ D 
from (||). The contribution of these terms, with the phys- 
ical cutoff at lower distances, would produce a smooth 
function added to (|l2|). This cutoff-dependent contribu- 
tion, being proportional to Np and finite at q — 0, may 
be ignored here. Note that this procedure of subtraction 
does not work when II(r) is a polynomial in 1/r; that is 
why we are left with the divergent expression for II (q) in 
even dimensions. We discuss the physical case D = 2 in 
more detail below. 

For the dimensions D = 1 and D = 3 the Eq. (O) reads 



as 



U(q) = N F tanh -1 \/l - k 2 



1-k 2 



(13) 



and 



U{q) = N F y\ - k 2 tanh" 1 yjl - k 2 - 1], (14) 
respectively. We plot these functions in the Fig. [l| 
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FIG. 1. The q— dependence of the Cooper loop. 

For completeness, we consider now the case of empty 
conduction band, i.e. [i < 0. It is convenient to use here 
the notions similar to the above Fermi momentum and 
the density of states : 



hp 



2m\fi\, N F 



2ttTW- 




A calculation similar to the above one gives : 
n(r) = \n\N F T 2 [v + l]V(k F r) 



= 2 [2/x 



(2a:), 



(15a) 
(15b) 



which shows the exponential decay of correlations at R 3> 
kp 1 . The counterpart of ( p"5| ) in the (/—representation 
(with the same restrictions) is found in the form : 



Il(q = 2k F n) = - 



Np 7T 

2 sin ttv 



[1 



(16) 



Again, one sees from (|I5[), (16) that in even dimensions, 
smitv = 0, the expression ( |16| ) is infinite. 

Now let us closer analyze Tl(q) in two dimensions and 
return back to (|l(]). In the realistic situation, the range 
of validity of the ellipsoidal form of dispersion d3|) is re- 
stricted to: the inverse interatomic distances, a~ . It fol- 
lows thenO that the range of applicability of (||) and, con- 
sequently, of (£t) is restricted from below by r > a. At 
larger scales the decay of H(r) becomes exponential one 
due to finite temperature or lifetime effects. In view of 
the mean-field equation for the transition temperature, 
cited above, we discuss the former source of exponential 
decrease of n(r). 

The main difference of the case of finite T is the sum- 
mation over the odd Matsubara frequencies lui in (|^), 
instead of Gaussian integration. This summation yields 
the Jacobi 9— function and we are unable to proceed 
further towards some analogue of (||). It is clear how- 
ever that the modification of (||) at large distances be- 
gins when one is no longer capable to ignore the discrete 



character of the Matsubara sum. In its turn, this is antic- 
ipated when the smallest Matsubara frequency provides 
the exponential smallness of the Green's function (0). 
The latter condition takes place at the distances r such 
that Re(y/-2mr 2 (fi + inT)) > 1. The simplest way to 
see it is to use the integral representation of the func- 
tion Kq 2mr 2 y//i 2 + uj'f jt in (||) and to integrate over 
t, then one obtains 



II(r, T) = 2TN F 



dOK y/h cosh(9-fo , (17) 



with b — Aymr 2 , b\ = Amr 2 \J + <-o 2 . The main contri- 
bution to (^) is given by 9 ~ 0, when the argument of 
Macdonald function \Jbi — b — 2Re{y— 2mr 2 ((i + iu>i)). 
A simple analysis will convince the reader that at large 
distances (i.e. when y/bi — b > 1 for all I) one can leave 
only the first terms u>i = ±nT in the sum (p"7j). The 
asymptotic behavior is given by (cf. (|To|) ) : 

n(r,T) = ^e-^0(l), r>£, (18) 

where the temperature coherence length £r depends on 
the ratio T//z and could be written as (cf. Ref.t3) 
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mm 



(19) 



Thus the conventional scale £t oc T^ 1 is changed by the 
shortet-one £x ck T -1 / 2 when the chemical potential is 
small.Eil As a result, the Fourier transform II(q, T) is al- 
ways finite and 



m 



n(? = o,r) = ^io g 



-o(i) 

a 



(20) 



In contrast to the three-dimensional case, Np is con- 
stant in 2D and H(q = 0,T) remains finite as kp — > 0. 
In this limit the dependence £t ~ T -1 / 2 only halves the 
usual value of large logarithm. It means that in two di- 
mensions the superconducting instability is possible at 
the arbitrary small filling ratios. Given a proper sign of 
interaction, even two electrons can form a bound state. 
This statement is distinct from the BCS picture and cor- 
responds to the fact that the bound state in two dimen- 
sions is possible even in the infinitely shallow quantum 
well. 

Let us discuss the possible generalizations of our main 
result (||). It would be tempting to obtain a closed ana- 
lytical form of the dynamical Cooper loop II(q, Q). Some 
analysis shows however that it is difficult, if possible. 
The main obstacle is the inapplicability of Eq. (||) , since 
for all Q there exists u> such that the requirement 
Arg(a + b) < ir/4 is violated. On the other hand, the ob- 
tained formulas could be generalized towards more real- 
istic anisotropic dispersions, in the spirit of consideration 
in RefJl 
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Concerning the possible implications of the SC chan- 
nel anomaly in two dimensions we wish to point out the 
following. In one of the possible scenarios of the physics 
of the high-T c cuprates, one discusses the formation of 
the spin-density wave state in these systems.!!! As a re- 
sult, the Luttinger theorem is violated and the Fermi 
surface can consist of small (ellipsoidal) pockets centered 
at (±7r/2, ±7r/2) points in the Brillouin zone. Our first 
observation here is that the small volume of these pockets 
does not influence the tendency to SC pairing. Second, 
the above factor 1/2 before the large logarithm in 11(g) 
is easily compensated by the number of pockets (four in 
the considered case). Next, one can see that the inter- 
ference between the electrons belonging to the vicinities 
of ("7r/2, ±7r/2) and (— tt/2, =F7t/2) ( "intervalley" scatter- 
ing) gives the same logarithmic enhancement for Tl(q) at 
q = (tt, 7t). This last possibility could be of special in- 
terest, giving rise to the non-local SC order parameter 
(d— wave pairing) and to the interference with the mag- 
netic AF channel. The consideration of these effects is 
however beyond the scope of this study. 

Summarizing we found the r— dependence of the static 
Cooper loop in any dimensionality. Its counterpart in 
q— space has the logarithmic singularity at q = and no 
singularity at q = 2kp in any D. In two spatial dimen- 
sions the bound state of two electrons is possible even 
in the limit of the empty conduction band. Although 
this fact is enough well-known, it should be taken into 
account in constructing the theory of high-T c cuprates. 
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